General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



f 

\ 


i 


f 


l^IVEUSHT OF CALIFORNIA. Rl-'JIKELEY 


hkhi:ei.ky • PAVis • luvisi: • i.os angeles • hiveksidk • san piu;o • san utANCisro \ 



E7 6- 1 0-3 7 


SANTA IIAHIIAliA • SANTA CWUZ 


NASA CR- 


SPACE SCIENCES LARORATORY 


BEREEEEY, CALIFORNIA $>4720 






TWX: UC SPACE BERK 
(910) 366-7945 

Remote Sensing Research Program September 17, 1975 

260 Space Sciences Laboratory 
University of California 
Berkeley, California 94720 


>r 

'M 

C 

in n 

'SI O' 

I _i ^ 

o u m 

r- z a 

Z 3 0 


-n 

'V 

'J5 


'NJ 

O 


> T ~ _l 

'J • O 

>-y. < Q. > vO 

Q. UJ JJ O 

5 X at i: 

— i o 

oO 

u. 3 n < 

LJ H- kJj — 

►— o C a. 

Z 'U 2 3 



D Z >- rx 

-J 3 -I -) 

O 'L — 

> OJ i- 

Li. ►- t— n 

I AJ < w" 

> > 3 -< 

3 -H O 

3 *— • 

»- *1 2 
t/1 ►- 2 ' 3 
-*J 3 -4 O 

^ a -H • 
ac <3 -< 't 
lU O “i 

t t 

0 3 i- _/ 

1 ~ t- ' 3 

o X a. 

jj 

i- i: -a 3' o 

w a. — — < 7 


Dr. A. E. Potter, Technical Monitor 
NASA-Johnson Spacecraft Center 
Mail Code TF 3 
Houston, Texas 77058 

SUBJECT: Quarterly Progress Report, Contract 

Contract for: Study-Development of 

techniques to IVheat Identification. 

Dear Dr. Potter: 

The objectives, approaches, supporting subtasks, and milestone schedules 
for the U.C. Berkeley LACIE tasks were given in the January 1975 proposal 
to JSC and refined at the beginning of the contract period (in conference 
with J. Deitrich and J. Garcia) as documented by the copy of the enclosed 
transmission (enclosure #1) to J. Garcia dated June 18, 1975. 

The major activities engaged in during the first 90-day performance 
period (May 15, 1975 to August 15, 1975) for each of the three tasks within 
the contract were as follows: 

TASK I ( 4.1.1.1a) Evaluation and Revision of Techniques for Manual Inter- 
pretation 

A prime objective of this task is to quantitatively evaluate various 
LANDSAT image sets in order to determine the optimum sequential data set(s) 
needed for accurate identification of the wheat crop. In addition, this 
task will determine expected accuracies for interpreters using optimum and 
selected less than optimum (missing data effect or desire for early season 
estimation) sequential LANDSAT data sets. The activities of this first 
reporting period centered on data acquisition and preparation, and area 
familiarization. The specifics of these activities are described below. 

Data files for LANDSAT imagery, ground data, and ancillary data were 
set up for the three agreed upon study regions - Kansas, North Dakota, and 
Montana, as well as for other possible regions of interest. 
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All data received from the J.S.C. technical monitor as well as that collected 
through our own efforts has been catalogued and filed for easy assess and to 
insure against loss. 

A library search for literature pertaining to the geography, agricultural 
environment and agricultural management practices of the study regions was 
undertaken to provide ancillary data for study site familiarization of persons 
who would be involved within the overall project and specifically to provide 
adequate background material for training of interpreters involved in the 
testing phase of this task. The bulk of this activity has been completed for 
Kansas and will proceed for the other two study regions (North Dakota and 
Montana) in a phased manner relevant to the milestone plan for Task I. 

Additional ancillary data in the form of weekly, monthly and annual 
reports from the respective states' Departments of Agriculture in cooperation 
with the U.S.D.A. Statistical Reporting Service have been subscribed to and 
are currently being received for all study regions. 

Based on initial work on the five intensive test sites in Kansas, the 
data found to be most useful as background material and as area familiarization 
material have been: 

1. The sequential full frame LANDSAT I color infrared color 
composites of the test site regions; 

2. The U.S.D.A., Soil Conservation Service Soil Surveys for the 
counties within the two crop reporting districts considered 
(the southwest and the central crop reporting districts) — 
particularly the sections of the reports which discuss the 
general agriculture of the county and the use and management 
of various land use classes such as dryland areas, irrigated 
areas and rangeland areas. 

3. The weekly weather and crop reports for Kansas, the monthly 
crop production reports for Kansas, the annual crop acreage 
summaries for Kansas 1971-1974, and special reports on grain 
quality and grain varieties for Kansas published by the 
Department of Agriculture and the U.S.D.A. Statistical Reporting 
Service. 

The five Kansas Intensive Test Site areas have been enlarged and copied 
onto color negative film for all dates of LANDSAT I color infrared composites 
received to date. For eight dates of LANDSAT I coverage of the Morton County, 
Kansas ITS for which we have not received color composites, but for which we 
have received the separate black and white positive transparencies, color 
infrared enhancements were made and enlarged to complete the sequential color 
infrared coverage for this site. The enlargements of the five Intensive Test 
Sites will be presented to the interpreters as part of the testing phase of 
Task I. 
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Analysis of the ground data, ancilliary data, and sequential imagery 
for the southwest and central crop reporting districts of Kansas is 
proceeding so that some preliminary screening can be done before the 
final design of the interpretation tests, and so that an efficient 
training package can be presented to the interpreters. 

TASK II (4.1.1.2f(l) and 4.1.1.2f(3) Development of a Spectral Signature 

Extension Methodology 

The primary objective of this task is to develop wheat analogue area 
stratification strategies and spectral signature prediction models to 
enable accurate, cost-effective wheat acreage determination in LACIE sample 
segments lacking ground data. 

As specified in the contract, there are five major supporting subtasks 
to achieve this objective. Progress on each task for this progress report 
period is cited below. The overall U.C. Berkeley approach to the signature 
extension problem is summaried in enclosure no. 2. An updated supporting 
data request for this task, dated July 8, 197.S and sent to J. Garcia, is 
included here as enclosure no. 3. 

SUBTASK I: Spectral Signature Extension Methodology Design 

The first component of this subtask is the specification of a stratifica- 
tion strategy to delineate wheat analogue areas. A wide review of the 
literature (including the following important references: Pascale and 

Damario 1902, Peterson 1965, Klages 1942, and Nuttonson 1956) combined with 
U.C. Berkeley's previous work indicates that the following variables should 
be used for stratification: 

. General soil type 

. Land use/cropping practice 

. Average growing season day-degree sums (April -June) 

. Avertage growing season precipitation (April-June) 

. Average last date of spring frost 

. Average temperature and/or average minimum temperature for the 
coldest month of the year 

These stratification variables are then used for delineation of analogue 
areas in one of four "first-case" ways that have been developed during the 
reporting period. Basically, these stratification techniques involve the 
following procedure. First, a single time of year is identified when soils/ 
land use/cropping practice patterns of interest are most separable on LANDSAT 
false color 9x9 inch format transparencies. For the crop reporting district 
or state of interest, these LANDSAT transparencies from the specified period 
from a given year are then enlarged (approximately 1.5 times) to print form 
and mosaiced together. 
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Clear acetate is placed over the image mosaic. County and other 
crop reporting unit boundaries are then drawn on the acetate relative 
to the image according to available planimetric map references. With 
the aid of general soil association information (e.g. regional soil maps 
or county SCS soil reports) general land use/cropping practice descrip- 
tions and/or two dates of supplemental LANDSAT transparencies that bring 
out different aspects of soil- land use characteristics, a photo-interpreter 
delineates on the acetate overlay the areal extent of the soil/land use/ 
cropping practice patterns as seen on the LANDSAT imagery. 

Onto the soil/land use/cropping practice mosaic is overlain broad 
threshold lines of the temperature and precipitation stratification variables 
cited earlier. Generally, the isolines of the temperature variables are 
very positively correlated with one another, so that only one or two need 
be considered for stratification. 

The resulting soil/land use/temperature/precipitation stratification 
defines unique wheat analogue area en.:ironments . Each of these environments 
should give rise to similar sequences of wheat spectral signatures during 
the life cycle of that crop, no matter where that analogue area type occurs 
in the world. Hence, variance among wheat signature and confusion crop/wheat 
signature mixes should be minimized , thus tending to maximize the success of 
signature extension models employed within each analogue area type. 

The general signature extension (prediction) model selected for experimen- 
tation is a regression-like equation known as an analysis of covariance model 
employing both qualitative (e.g. crop type) and quantitative (e.g. a continuous 
variable such as temperature) variables. A regression equation will be employed 
as the "first case" expression of the more general model. 

The regression equation will have as its objective the prediction of 
the mean spectral response for a given LANDSAT band for a given crop residing 
in a recognition segment. The "prediction" or "independent" variables in the 
regression equation used to predict the dependent variable, (mean spectral 
response) have been selected to be: 

. Crop maturity codes from A.I.s 
. Soil color codes from A.I.s 
. Growing season day degree sums 
. Planting season day degree sums 
. Growing season precipitation 
. Planting season precipitation 
. Sun angle 
. Scan angle 

. In-segment LANDSAT spectral data such as the BAND 7/BAND 5 ratio 
(crop calendar correlated) 

. Atmospheric transmittance (generated from ground horizontal visibility 
data obtained at airports) 

The regression coefficients are estimated for these predictor variables from 

data in training segments for wheat and major confusion crops for each phonological 

phase. 
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Signature predictions can then be made for LANOSAT bands 4, 5, 6, 7 for 
a given phonological phase on a field basis or for clusters identified 
within the recognition segment through non-supervised classifiers. Crop 
identification is then attached to the fields and/or clusters by statis- 
tically testing for the similarity between the predicted signatures (assuming 
a given crop - one regression equation set per crop per phonological phase) 
and the actual signatures in each LANUSAT band for the given field or cluster 
in the recognition segment. Statistical tests of similarity include single 
band students-t tests and multiple b=»nd multiple comparison tests such 
as those given by Scheffe (1959) . 

The crop type having the least significant statistical difference 
between predicted versus actual spectral signatures is assigned to the 
field or cluster in question in the recognition segment according to mean 
and covariance statistics for crop types obtained from the individual fields 
or in the statistical clusters in the recognition segment identified by 
crop type in the previous step. 

SUBTASK 2: Acquisition and Preprocessing of LANDSAT and Supporting Data 

Work to date has focused on the southwest and central crop reporting 
districts in Kansas. All LANDSAT 9x9 inch transparency data have been 
obtained for these areas for the 1973-74 crop year except for the color 
infrared composites of Morton Co., Kansas as noted in the Section on Task I. 
1974-75 crop year LANDSAT transparency data acquisition is as yet incomplete. 

JSC digital data for the LACIE 1973-1974 Kansas R5D segment set is on 
hand with the associated crop I.D. and ancilliary information. Similar 
digital SRS segment data and associated descriptive information for the 
Kansas 1973-74 crop year is now arriving from JSC. U.C. Berkeley is presently 
reserving a request for full or partial scene LANDSAT digital data for a 
portion of Kansas to allow a broader segment population on which to test 
its technique. 

Regional soil map data for Kansas has been obtained and the published 
county-wide soil survey data (USDA S.C.S.) for the two crop reporting districts 
is on file. Climatic ground station data for the 1972-1975 time period 
for Kansas has been obtained from the U.S. Department of Commerce, Environ- 
mental Data Service. Airport ground horizontal visibility data for select' 
dates 1973-1974 in Kansas is now available. Historical data, crop calendar, 
and cultural practice information to be used in TASK II parallels that 
obtained and described for TASK I. 

Data for the complete western two-thirds of Kansas desired for stratific- 
ation by our LACIE technical monitors is incomplete with respect to LANDSAT 
'73- '74, '74- '75 transparency coverage. Most of this data should be acquired 
for stratification purposes in the next two months. 
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Most data for the North Dakota Study Area as well as for the otlier to 
be defined study area (another state or portion of a foreign country) 
is not yet available. The exceptions to this statement include complete- 
ground station climatological data for the period 1972-1975 for North 
Dakota, Texas, Illinois, Indiana, Montana, Washington, and the world as 
published by the Department of Commerce. Also significant historical, 
land use, cropping practice, and general soil/climate data is presently 
being obtained for the Soviet Union. The Russian data is being obtained 
for a tentative U.C. Berkeley LACIE TASK III effort, but could be used 
for TASK II signature extension research as well. 

SUBTASK 3: Stratification for Analogue Areas 

Stratification using one form of "first case" stratification strategy 
discussed previously will soon (2 weeks) be completed for the southwest 
crop reporting district in Kansas. Completion of first-case stratification 
of wheat analogue areas should be completed in the central Kansas crop 
reporting district soon thereafter. Stratification variables receiving 
the most significant attention at this iteration are general soil/land use/ 
crop practice strata, growing season (first-case definition -- the three 
month period following the month of the vernal equinox, i.e. for Kansas, 
the April to June period), day-degree sum, and the growing season average 
precipitation. The average last-date-of-frceze and the vernalization 
factor of average temperature or average minimum temperature in the coldest 
month of the year (correlated with indices of number of days below 5 C., 

etc.) are being considered for inclusion in the stratification scheme. 

Stratification for the remainder of the western two-thirds of Kansas 
has been limited primarily due to the lack of several LANDSAT 1 transparencies 
which have been requested through our TASK II monitor and are, as we under- 
stand, forth-coming, and to manpower limitations. 

Stratification for wheat analogue area definition in other states (countries) 
has not yet begun as per previous milestone schedule guidelines. 

SUBTASK 4; Determination of Numerical Constants in the Spectral Signature 
Prediction Equations 

Predictor variable values arc being calculated for 1973- '74 LACIE sample 
segments in Kansas. For a simulated set of training segments, regression 
coefficients for each prediction variable will be estimated according to 
least squares analysis of paired LANDSAT band, mean spectral data and corres- 
ponding values of signature prediction variables. 

Regression constants will be estimated for each regression equation 
needed to represent each important combination of crop type and phcnological 
phase. These regression equations will then be applied to fields and clusters 
not used to develop the regression coefficients, but for which the true 
crop identity is known, to confirm that each equation is accurately predicting 
the expected signature for the crop in question. Wlien predicted and actual 
signatures for the same crop do not match according to the statistical 
similarity tests, then the training segment and training set selection procedures 
will be modified. 
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Regression coefficients are being determined for the southwest and 
central Kansas crop reporting districts as first priority. These two 
crop reporting districts contain all five intensive test sites in Kansas. 
Temperature, precipitation, and transmittance prediction variables are 
being calculated for a sample segment as a whole according to a weighted 
average of values from the nearest ground meteorological stations. The 
weight for a given station value is inversely proportional to its distance 
from the segment center. Other data such as in-segment LANDSAT band ratios 
and crop maturity is presently being considered on a field or cluster 
specific basis. 

SUBTASK 5: Application of the Spectral Signature Prediction Equations 

to Corresponding Simulated or Actual Remote Analogue Areas. 

This last step in approach to signature extension awaits further 
development in the foregoing steps. As a first priority, signature 
extension by use of regression prediction equations within the same wheat 
analogue area type will be performed in and between the southwest and 
central Kansas crop reporting districts. Signature extension over larger 
areas of Kansas and other states will follow according to the milestone schedule. 

TASK III: (4.1.3.1b) Variable Probability Sampling for Acreage Estimation 

The objective of this Task is to develop and test photo interpretation 
techniques for definition of wheat sampling strategy. Essentially, this 
task involves the determination of the photo interpreter's ability to 
provide a quick estimate of wheat acreage in the population of sample segments 
in order to provide better wheat density information on v/hich to allocate 
sample segments for computer information. It is important to note that in 
this first sample stage photo-interpreter information can also provide a much 
more precise estimate of wheat acreage by reporting unit and potentially, 
significantly reduce the number of the much more expensive computer analyzed 
segments needed at the second stage (phase) of the sample. In other words, 
if quick-look photo interpreter estimates of wheat acreage by sample segment 
are significantly correlated with computer results and/or ground truth (R§D) , 
then significant improvements in wheat acreage estimation accuracy, precision, 
aid cost may be obtainable. 

Work to date on this task has concentrated on selecting areas for testing 
and demonstration. Two crop reporting districts in Kansas, the southwest and 
central, have been sleected for the winter wheat type. Sample grids (a matrix 
of 5 X 6 mi sample segments) are being constructed for overlay onto LVNHSAT 
transparencies and signature extension LANDSAT print mosaics. The sample grids 
themselves are being constructed for the state of Kansas. The portions over- 
laying the two crop reporting districts mentioned will then be used for plioto- 
interpretation . 

Photo interpretation will proceed in the first-case demonstration on early 
July 1973 LANDSAT imagery. This single growing season date for western 
and central Kansas appears to contain a maximum of wheat acreage information 
with a minimum of crop confusion information for the quick look approach. 

At this time of year in this winter wheat area most dryland and irrigated 
wheat fields have a distinctive yellow to beige tone on the false color FJVNDSAT 
data. Other regions may require more tluin one date of 1.ANDSAT imagery for 
optimal results with the "quick-look" aporoach. Furthci- research on the Kansas 
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areas could also indicate a multidate requirement. 

Photo interpreter training Ueys/aids are presently being developed 
for the two Kansas crop reporting districts ci ted earl icr. These keys 
are based on cropping practice/soil relationships and crop calendar 
information as they relate to the image features of wheat on the LANDSAT 
data selected for interpretation. The training aid information is being 
derived from data sources cited in U.C. Berkeley Tasks I S II. 

Demonstration sites for areas other than Kansas are yet to bo finally 
determined. Recommendations from our technical monitors indicate that 
Montana or Idaho would be good regions in which to demonstrate the "quick- 
look” wheat acreage estimation technique for a mountainous area where small 
agricultural areas are found within interfingering valleys. But the decision 
as to the final area of study within one of these states awaits an evaluation 
of the availability of cloud-free properly timed with crop phase imagery. 
Evaluation of available imagery is now in progress. 

Foreign country discussions have centered on the Soviet Union, partic- 
ularly one oblast or a portion of an oblast. U.C. Berkeley has been studying 
Soviet Union historical wheat data, climate, soils, and existing LANDSAT 
dr.ta to determine the location and relative concentrations of winter, spring 
and mixed wheat in that country. The purpose of this Soviet analysis, now 
complete, was to enable U.C. Berkeley to make a recommendation to our LACIF. 
Task III monitor concerning the location of "quick-look" wheat-acreage demon- 
stration sites. 

The recommended Russian areas are either a mixed wheat area in 
the region defined by the cities of Penza and Kuyushev (north and NNIV of 
the Caspian Sea) or a spring wheat area in northern Kazakh S.S.R. NE of 
the Caspain Sea). The former area is an established small grains producing 
region while the latter is located in a "new" agricultural lands area 
brought into production by irrigation. The latter region has a climatic 
regime apparently similar in many respects to the Southwest U.S. Variability 
in wheat acreage from year to year in the new lands area may be quite high. 
The decision as to which region should be used for demonstration of the 
Task III "quick-look" acreage estimation technique depends on tlie sensitivity 
of the LACIO acreage estimation formulas to the two types of Soviet wheat 
producing situation cited. 

These foreign country recommendation study results, just obtained, have 
yet to be discussed with our Task III technical monitor. 
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If there are any questions concerning this report, please feel free 
to give us a call. 


Sincerely yours, 



RSRP-LACIE Project Leader 
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Randall W. Thomas 
Technical Services 
Branch Manager 


cc: Messrs. John Deitrich 

Joe Garcia 
A1 Feiveso 
Joe Downs 
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Hvaluation and Itcvision of Techniques for Manual Interpretation 
1.1.1 Objectives 

The primary objective of tliis ^)roposcd task is to evaluate various 
manual interpretation techniques for the identification of wheat on a world 
wide basis and to revise them as necessary. Specifically, the techniques 
will entail the human analysis of L/\NDSAT iinajicry and ancillary data for the 
purpose of providing reliable training data for nutoinatic classification 
systems in remote areas where ground data arc not available. 

The emphasis of the proposed manual interpretation work will include 

the following subtasks: (a) the quantitative evaluation of various I.ANDSAT 

# 

image sets, and (b) the definition of the best procedure for tlic image analyst 
to use to identify wheat. Note that these two subtasks arc closely related. 

For example, on optimum set of imagery to be tested in Subtask A nay have been 
defined on the basis of temporal information determined in Subtask B. 

Initially, however, these two subtasks, will be evaluated as independently as 
possible of each other. 

A. Quantitative Evaluation of lANDSAT Iiiiai jcry 

Because of the numerous multi temporal image sets and riultiband/ 
multidate image combinations that can be produced from LA'vUSAT data, llic 
determination of the best set of temporal images needed fur accurate crop 

identification will be a major objective of the proposed study. 

. r 

Specifically, the interpretation accuracy tliat can be cj^pocLcd fio;:! inter- 
preters using 1) optimum photogenic dale or combination of da.tcs (sequential 
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technique) durinc tlic Rrowinp season; nnd 2) less than opt inuinrplioloi^cnic 
dales or combination of dates durjnt; the p.rowin}* season (i.c. other than an 
optiiiuini date due to mission images as a result of cloud cover, etc. or to 
incomplete optimum image set in an attempt to make estimates of wheal acreage 

early in the growing season) will be determiried from interpretation tests. 

B . Image A nal ysis Tecbnicpics 

The objective of tliis sublask is to define and evaluate those 
materials, equipment, and techniques which arc required by the image analyst 
to complete the task of identifying wheat. This will include, in part, 
determining the utility and optimum mix of the following materials and 
techniques; 

1. reference materials -- crop keys, crop calendars, data f)’om 

/ analogous areas, etc. 

2. auxilliary data — high altitude photos, low altitude photos, 
historical agricultural data, spectral data, meteorological 
data, etc. 

3. convergence of evidence -- field size, cropping pattern:-, cultural 
features, etc. 

4. conference system -- joint interpretation with two or more 
analysts. 

1.1.2 Approach 

Micrcvcr possible, a quantitative cvahuilion of all the sublasks will 
he wade, in order to determine how the accuracy levels and associaloil error;, 
vary with tl»c different data sets that arc examiiud. It must be cmplia sized 
that in order for any of the results to have meaning witli the LVClIi, all 
existing ground data, image ;;ets, crop calenders, cropping practice infoim.-tt ion, 
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anti liANDSAT computer compatible tapcE, ir>t be made available to the Ibanote 
Sojsinu Research Rrocram, Not only will this make the results more moaninyful 
with rcfianl to LACIH, but it will minimize duplication of ciTort between 
participatinj; organi zations. 

A . Quantitative nvalu.iti on o f l.ANDSAT J ma;>cry 

Image intcrj)retation tests will be given to determine I lie exjiectcd 
interpretation accuracy of optiinur.i image sets for the identification of v.iicat. 
These tests will be given to skilled image analysts from both the Rfild’ and 
the LACin Lockheed support gr 9 up. Test results will be analyzed using an 
analysis of variance design so tliat tlic levels of accuracy attainr.Iilc from 
the different sets and their conficloncc limits can bo dcleriiLiiod. 

Most of the image sets that will be tested will bj produced optically, 
photogrifphically , and/or clcctroni cal ly by personnel fioi:i the .lolinson 
Spacecraft Center and the RSUP. In order to completely evaluate the state-of- 
the-art in image production, however, a limited number of image sets will be 
purchased from educational facilities and co:mncrcial organ Lzntiofis for testing 
purposes. 

B. Image A nal ysis Techniques 

Hie results of tlic image interpretation tests given to evaluate tl’.e 
image sets will bo used also to evaluate the image analysis process, in order 
to determine Cl) the contributions to identification accuracy made by 
various crop keys and ancillary data, and (2) the effect that nirsing data 
have upon usefulness of these data sets. 

1.1.3 expe cted Resu lts 

Tlio cxjiccled results include the document at i on , evaluation, and revi;.ion 

of tochn'Kine.. '’'’r ilie manual interpretat itm of LAN’Di’iAT s;.aimev il.tin foi* the 

purposes of training a comjnjt or ;i 1 j-ori t Iim to ideiiti fy i.dicat. 
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llicsc rosullr. v/ill incliulo; (1) a ilc:.Lri pi ' on o(* the opt it.i in LAMDSAT 
sets for ident j fic.'itJoii of whcrit alonp, v/ith tlio accuracies of 
interpretation expected usinj; the optiinnni sets as i;ell as less than opt ii.iiaii 
sets, and (2) a description of the necessary and r.uffi c i cut ancillary data 
needed by the interpreters. 

1.1.4 Mile ston e Tla n 

The inilostono plan for Task I is p.ivcn in Table 1, It sh.oulil be noted 
that all time lines vdll shift to the rij;iit pro]iortional to the delay in 
receiving bANIlSAT and auxilliary data. 
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TARLU 1. Mll.l'STONU PLAN 1-OR L\CIl; TASK I 


RflPBODUCIBtLrrY OF THE 
ORIGINAL PAGE B POOR 


Subtask A, Qiirmtito t ive nv.-iluation of Lai idsat Tmn.f»ory Sot(s). 


1. Dctcrim' nation of optiiiium sequential data sct(s) for Kansas test 
sets (winter wheat areas). (Data expected from JSC by 7/15/75.) 


a. Rank all dates according to their effectiveness in identifying 

wheat. (10/15/75) 

b. Description of optiinua sequential data set(s) for Kansas test 

sites - best 4 or 5. (10/15/75) 


c. Interpretation accuracy expected from optiinu-! data set for 
Kansas test site as determined from tlic interpretation tests. 

(10/15/75) 


d. Interpretation accuracy expected from less tluin optimum data 
sets for Kansas test sites. (10/15/75) 

2. Determination of optimum sequential data sct(s) for North Dakota 

test sites (spring wheat areas). (Data expected from JSC by S/1/75.) 


a. Rank all dates according to their effectiveness in identifying 
wheat. (12/31/75) 


b. Description of the optimum sequential data set(s) for North Dakota 

spring wheat sites. (12/31/75) 

c. Interpretation accuracy expected from optimum data set for 

North Dakota spring v.-licat test sites. (12/31/75) 

d. Interpretation accuracy expected from less tlian optiir.iui sequential 
data sets for Nortli Dakota spring wheat test sites. 

(12/31/75) 

3. Determination of optimum sc(iuontir.l data sct(s) for an area of 
mixed winter and spring wheat production (possibly Montana or 
Soutli Dakota). (Data expected from JSC by 8/1/75.) 


a. 


Rank all dates according to their effectiveness 
winter and s])rir.g wheat. 


in identifying 
(5/1/76) 


b. 


Description of oplimun sequential data sct(s) for an area of 
mixed spring and winter wheat production. (3/l/7i)) 


c. Intovprtaal ion accuracy expected from optimum sequential data 
set for an area of mixed sprinj; and winter wheat j)ro(liict ion 

(3/1/76) 

d. Interpretation accuracies expected from less tluin optimum data sets 
for an area of mixed spring and winter wheat production. 

(3/1/76) 

Suh I a s k B. Evaluation o f Image Analysis Techniques 

Description of necessary and sufficient data nputs (imagery and ancil- 
lary data) for the interpretation of winter wheat based on the Kansas data, 
spring v/heat based on the North Dakota data, and mixed spring and winter 
v;heat based on the data of an area of mixed product ion- possibly Montana or 
South Dakota. (5/11/76) 

a. Large area full frame input to the interpreter vs. small area 
test site only (5 X 6 mile segment) input to the interpreter. 

(5/M/76) 
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1.2 Task II: Development of a Spectral Sif.natnrc Extension Mothodolo|;y 

1.2.1 Objectives 

The primary objective of this task is to develop analogue stratification 
strategies and spectral signature prediction models tliat v/ill alio;-/ accurate, 
cost-cffcctivc classification of v.’heat areas in agricultural regions of the • 
world lacking current ground data. 

The development of this spectral signature extension methodology will 
involve five supporting subtasks: 

A . S pectra 1 Signature F.xtension Mot liod ology De sign 

This subtask will involve (1) development of stratification 
strategy for analogue wheat areas, (2) selection of the spectral signatui-c 
prediction model to be employee! and a definition of variables utilized in the 
model, and (3) design of supporting sampling strategies. These design specifi- 
cations v/ill guide and coordinate progress on the remaining subtasks. 

B. Acquisition and preprocessing of LAN DS AT a nd Supporting Data 

This subtask will include acquisition and indexing of LA'>!DSAT data 
and supporting ground and aerial data, as well as meteorological, soil, and 
plant culture information as required by design specifications for stratification 
and model variables. Transformation of data planes to a conunon kA?.'i)SAT 
coordinate system and sul)sequcnt information merging will also be performed 
where necessary. Integrated information management, as this stage will 
permit most cfficioit processing in the remaining thi'cc subtasks. 

C. Stratif ication for Aivi lo guc Arens 

Application of the criteria for delineating analogue areas over 
specified regions will be performed with subsequent refinement of stratification 
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mcthodolofiy. AnaJoijiic areas arc defined as ajiriculLiiral lands potentially 
giving rise to similar wheat signatures (spectral and spatial) over the 
life cycle of wheat due to similar climate and soils conditions as well as 
similar cultural practices. After strafification, each agricultural stratum 
will define a particular wheat analogue area type. Within a given wheat 
analogue area type there will exist non-remote analogue areas (current ground 
data available) and remote analogue areas (no significant current ground 
data available). Stratification of regions for analogue areas is designed to 
significantly reduce variation in predicted spectral signatures and also the 
cost of signature extention and stratification. 

D. Determi nat ion of Numerical Constants in the Spec tral Sig nature 
Prediction equations 

> For wheat and non-wheat training fields selected by image analysts 
in non-remote analogue areas parameters will be estimated for the spectral 
signature production equation for each analogue area type. Tests for statis- 
tical significance of variables included in each equation will then be conducted. 
A determination of the accuracy of the signature prediction models within 
non-remote analogue areas will allow further refinciacnt of model parameters. 

The objective of this subtask is to obtain the ability to generate accurate 

9 

channel means and variances for each wheat analogue area type under varying 
conditions. 

n . Applicat ion of the Spectral S i gnature Prcd i ctit)n l'(|iiations to 
Cor respon ding Simulated or Actual Remot e Anal ogue A reas 

Within remote analogue areas, data from wheat and non-wheat fields 
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located by imaj;c an Jysts will be entered into the ypectval r.in'witurc predic- 
tion equations particular to those strata. Predicted spectral estimates will 
then be statistically compared by channel to the corresponding actual 
remote area spectral signature means. If statisti cal ly significant differences 
occur in critical channels, then retraining for wheat and non-wheat locations 
in remote areas will be initiated. Otherwise, wheat areas will be classified in 
the given remote analogue areas according to spectral statistics particular 
to each remote area's training fields. Finally, confidence intervals around 
wheat acreage estimates in analogue areas will be constructed and cheeks on 
classification accuracies will be performed where possible. 

1.2.2 Approach 

The five subtasks outlined in Section 1.2.1 for the dcvclcpmont of 

> 

an accurate, cost-effective spectral signatui'o extension methodology will be 
conducted in a coordinated manner both stepwise and in some instances simul- 
taneously. The overall relationship of the five subtasks is shovni in Table 2. 
Methodology design and data acquisition must necessarily be the overlapping 
first steps followed by analogue area stratification, model parameter 
determination, and finally model application to simulated or actual remote ttaalogiic 
areas for wheat classification. On time completion of all tasks will be de- 
pendent in a sliding scale manner on tlie on-time RSRP receipt of funds and 
data from lAClF. management. Specific approaches for each of tlic five subtasks 
will be as follows: 
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A. Spectral Signature HxtcMision Motliodolo^’/ Do.sJf>n 

Stratification strategy for analogue ureas will involve 1) a specifica- 
tion of important stratification variables incliuling soil uiul meteorological 
conditions, cultural practices, wheat types and confusion crops grown, and 
2) selection of human and huinan-autoiantic tcchnicjucs to allow rapid delineation 
of strata according to stratification variables of interest. An important 
point to note is that a given wheat. analogue area type must occur within the 
same agricultural stratum. This means that a given stratur.i will be disjoint, 
that is, it will occur in several separate places over the globe. 

The spectral signature prediction model has been tentatively selected 
as the analysis of covariance model. This is a linear probabilistic relation- 
ship that includes qualitative factors and their interactions as in analysis 
of variance models while at the same time providing for quantitative variables. 
Hence the statistical significance of several spectral signature modifying 
factors in a p-way layout can be cost-effectively and realistically tested 
after correction for regression on a quantitative string of other variables 
also imp.acting on spectral signature. The dependent variable in the analysis of 
covariance model will be taken to represent the predicted mean value for a 
given spectral and/or spatial infoianation channel for a givci\ analogue area 
according to tlic values of tlic signature affecting vax'iablcs existing in the 
given analogue area at a given time. 

The selection of which qualitative and quantitative f.actors and their 
interactions to include in the analysis of covariance model will dopciul on 
factors identified in the litcx*aturc, in RSUP research, and in forums with 



other LACII: signature extension working group members as having significant 
effects on spectral signature. Moreover, any factor or variable inclutlcd in 
the model may be complex. That is, raw temperature or humidity values may 
not be used directly, but some physiologically related quantities generated 
by deterministic relationsliips could be employed in the model instead. 

Supporting sampling strategics utilizing LAMDSAT, aerial, and ground 
data will be used to calibrate signatures wlicrc ground data are obtainable to 
improve the precision and accuracies of tlic spectral signature extension 
equations. These sampling strategies will be designed to be cost-cfficicnt 
and workable in real-time v/heat acreage estimation situations. Moreover, 
estimation of variables employed in the signature prediction models will be 
based on sampling plans involving cost-cfficicnt combinations of data. 

r 

B. Acquisition (and Preprocessing) of L\i\’DSA7 and Sui i portiiig Data 

Tlie approach to an organized and readily manipulatablc data system will 
involve firstly an acquisition and indexing of LAND5AT digital tapes and 
transparencies for areas specified in the study. Digital tape reformating, 
calibration, study area extraction, and if necessary, generation of a conmion 
LANDSAT multidate coordinate system will follow. 

Acquisition and indexing of supporting aerial and ground data will 
he performed in non-remote analogue study areas and auxiliary information 
planes, including topographic, meteorological and cultural, will be obtained 
to tlic extent possible for .nil study areas. 

A generation of .nuxillary data planes providing information important 
in the stratification process and the analysis of covariance spectral 
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signature prediction model will then be performed and, when necessary, 
tronsforincd to the common I<ANl)SAT multidntc coordinate system. 

Finally the LANDSAT and auxiliary data planes will be merged where 
appropriate for stratification, model development, and wheat classification. 

C . Stratification of Analogue Areas 

Stratification for wheat analogue areas in specified study regions 
according to variables determined important in defining homogeneous wheat 
producing locations will be conducted according to manual or optionally manual- 
automatic technologies. Manual approaches will involve delineation of analogue 
areas on hardcopy false-color enhancements of IJVf^DSAT and superimposed 
auxiliary data combinations. Auxiliary data may include general climatic, 
soil fertility, drainage, and cultural practice indices. Human-automatic 
stratification techniques might involve human location of analogue area 
boundary coordinates on manipulatiblc TV monitor displays of merged LANDSAT 
climatic, soil, and cultural practice information. 

Analogue area stratification technology development for the period to 
December 15, 1975 will focus on the winter v;heat areas of the State of 
Kansas, particularly the designated LACIG intensive test sitc(s) v/hcrc 
special data for signature extension research will he collected by L/\CIG 
cooperatoi'S. Delay in RSRP receipt of data will cause delays in task completion. 

For the period December 16, 1975 to Hay 14, 1976 analogue area 
stratification tcclmology development will focus on spring wlicat j)roduction 
areas of North Dakota, particularly tlic designated L\Cin test 5itc(s) where 
Special data for signaturn extension research will be collected For LACIIi 
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cooperators. Again, delay in HSUI* receipt of data will potentially cause 
delays in task completion. 

D . Determination of Humev ical Co nst ants in the Spect ral Si g iiature 
Prediction Equations 

Image analyst selection of training fields for estimation of spectral 
signature prediction equation parameters (i.c. numerical constants) will be 
based on crop identification information gained in non- remote study areas. 

Klicre possible, crop identification information will be determined from 
analysis of LANDSAT, aerial, and ground data collected in cost-effective 
sampling schemes. Tlicse sampling plans may be the same as those designed for 
LACIE v/hcat acreage over non-remote analogue areas. 

Estimation of model parameters vdll proceed according to least- 
squares criteria for data derived from the geometrically coincident lANDSAT 
and auxiliary information planes occupying selected training field locations. 
Training fields will be located so as to provide information for a complete 
p-way layout with at least two samples per cell. 11ms all parameters for 
interactions of interest in the spccttal signature prediction model should be 
estimable. 

All parameters for the p qualitative (prcscncc-absonco) variables 
or their interactions in the analysis of covariance model will ho determined 
by calculating their values according to usual analysis of variance procedures 
and then subtracting from each such numerical constant a correction for 
regression on the quantitative • (continuous) variables not controlled in the 
p-way l.iyoiit. For instance, assume that the number of qualitative variables. 



p, is oqiinl to three, rurthcr assume that these three quali.ttrtivo variables 
consist of wheat aiialo|iuc area type, life cycle stage, and class of 
irrigation frequency. Then the effects (expressed as ninncrical constants) of 
those variables on wheat signature Avould be calculated according to usual 
analysis of variance rules and then corrected by regression of wheat signature 
on such quantitative variables as sun angle, inches of rain for a given time 
period, etc. In all cases the values of the dependent variable will be 
generated from the mean nuincrical values for given information channels on 
given training fields. In this way, the true importance or magnitude of 
contribution to a predicted spectral signature can he assessed for each combina- 
tion of qualitative factors after adjustment for the effect of raany freely 
varying environmental variables. 

t F-tests and multiple comparison tests including the Scheffe method 
will be made to determine the statistical significance of qualitative and 
quantitative variables and associated interactions included in the specti'al 
signature prediction equation for each wheat analogue area type. Those tenns 
not found to be significant will be dropped in order to strcoMline the 
prediction model. In this way a cost-effective model will be generated 
utilizing only data found to be significantly correlated with wheat signatures. 
Furthermore, estimates of variables important to signature prediction will be 
based on a cost-efficient sampling process utilized in the preprocessing subtask. 

Development of software for the analysis of covariance prediction 
model and associated tests of statistical significance for model parameters 
will bo designed for compntiI»ility with current data processing systems. 

A test will tlicn bo conducted in non-remote areas for the accuracy 
of wheat field signature idcntificatinn by tlio paiticular spectral signature 
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prediction model associated with each wheat analogue area type’. Where 
signatures for known wheat fields arc not accurately estimated, then further 
training field location and model parameter estimation will be conducted 
until statistically adequate signature predictions are made. 

E . Application of the Spectral Signature Pre d iction Equations to 
Corresponding Remote Analogue Areas 

Image analyst selection of training fields in remote or simulated 
remote analogue areas will be based on inferences of wlieat presence from 
various LANDSAT spectral and spatial enhancements and on auxiliary data such 
as crop calendars. Other LACIE work on remote identification of wheat will 
have important application here. Simulated remote areas will consist of non- 
remote locations treated as if they were remote. Aerial and ground di’ta for 
these simulated remote locations will allov; subsequent determinations of tlic 
accuracy of wheat classification after spectral signature extension. 

For wheat training fields selected in remote or simulated remote 
analogue areas, merged LANDSAT and au.\iliary data will be directly substituted 
into the spectral signature prediction equations particular to the given wheat 
analogue area types. The resulting signature estiraates by channel will tlicn be 
stati.stically compared through system compatible software with known corresponding 
values from remote areas. Tlicsc determinations of statistical significance 
will be tests of the null hypothesis that remote training fields arc in fact 
wheat. Candidate tests include siraplc pairwise t-tests between predicted 
versus knov;n channel weans and multiple comparison tests between several 
predicted versus knomi channel means. 
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• If statisticully sip.iiificant differences in key chmincls occur, 

then retraining for wlicat fields in remote areas must occur. This feedback 
to the remote wheat identification task will indicate whether further refine- 
ment of associated remote identification techniques is necessary for 
specific wheat analogue area types. Retraining will enable generation of a 
new set of wheat spectral signature predictions. 

If no critical statistically significant differences in channel 
means occur between kno\m and predicted wheat signatures tiicn wheat wi ll be 
classified in the remote areas according to usual methods. Specifically, 

t 

training statistics from fields identified as wheat by the spectral signature 
prediction equation will be used to classify wheat acreage in sample units 
for the particular remote analogue area. TTie bi^e, shape, and nu>iibcr of 

I 

sample uiiits selected per analogue area will of course depend on the wheat 
acreage estimation procedures used for I.ACIE. For most accurate wheat classifi 
cation, it may be necessary to identify non-wheat confusion crop training 
locations also by spectral signature prediction equations and then use the 
statistics for those non-wheat types in the classifier as other classes. 

Confidence intervals around wheat aci'cagc estimates will be 
formulated according to usual methods. It is important to recognize the 
overall significance of the confidence interval argument in the proposed 
spectral signature extension approach outlined here. Among other applications, 
confidence intcx'vals will be constructed ai'oiuid 1) estimated raw data variables 
2) estimated variables used in the spectral signature equations, ami 3) 
average channel values estimated by the prediction equations for remote areas, 
in order to determine the reliability of the ultimate acreage estimate. Tlicsc 
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confidence intervals v/111 also he used as criteria to adjust spectral 
signature prediction results when input data are lacking or inconplctc in 
remote areas. As such the intervals will constitute a portion of the 
sensitivity analysis for missing data. 

1.2.3 nxpcct ed Results 

The primary results from the proposed investignt ion v/ill include a 
documentation of techniques employed in the five suhtashr. and an cvaTiijtion 
of the timeliness, cost-cffcctivcness, and accuracy witli whicli wheat can 
be classified in simulated remote analogue areas. Any worknhlo t.f chnoJogy 
tliat is developed in tlicsc tests to employ spectral signature o:;l.cnsion 
methodology for large, rcal-tiiac global inventories should bo t fiinsfciab] c to 

user agencies. This technology ti. nsfcr will involve specif ic.;t i ons for 

r 

hardware modification in initial and more advanced spectral signature cxicaslon 
methodologies and will also include all software particular to Hie signature 
extension process. 


1.2.4 Milestone Plan 


• The milestone plan for Task II is given ir Table 3. 
that all time lines will shift proportional to the delay in 
and auxiliary data. •* . ». ;• , . 


It should bo noted 
rccc i V i ng I,A.‘!()SAT 
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TABLI! 3 MIM-STOfn I’lAN I’OU I.ACIII TASK II. 


1. Initinl SpectraJ Signature I-xionsion (SSP,) Melliodolnpy Di'sign. 

(Most of Mctlioilology Design be ilevolopcil iisiiij; Kansas «!aVa.) 

u. Development of a stratification strategy: variables, tedniiqiics, 

wiicat life cycle relations between nnilognc areas, and criteria 
for analogue areas, (nrief report or oral present alien for 
Kansas (9/l/7f.) 

final report due (irVlb/7I>) 


b. 


Selection of SSP prediction model 


lype(s) and variables 

the model - V/1/V.S 
tlio varialilcs -lO/l/yii) 


c. Design of statistical te.sts for varinbJe ineliniioa uilhin titc 
SSI: model (s). (5)/l!'./7r,) 


d. 


Specification of statistical decision algoritlv.i to determine 
significant differences betwcoji actual and predicted i j’.mi i ores . 

(0/15/7:0 


c. Design of SSH sensitivity analysis tests to (pi 'ni ify pr'rfoj-i.iatu'.o 
and isolate error sources. (I’rcl imiiuuy I’c.poj L • 10/15/710 
* I’inrjl report - 12/15/7.5) 

f. Design of SSH supporting sanpling strategies. (10/1/75) 

2. Acquisition of Laiidsat and supporting data for selected intensive 
test sites. 


a. Identify data retiuirements for Kansas (C/20/75) 

Identify data requirements for Mortli Dakota. (9/1/75) 

Identify data requirements for Other. (9/15/75) 

b. Acquisition of needed data for Kansas (7/15/75) 

Acqui.sition of needeil data for Nortli Dakota (10/1/75) 

Acquisition of needed data for Other. (10/15/75) 


3. Stratification of Analogue Ai'cas. 

Application of stratification criteria to selected I ■ ir 
a. Ktins.as (Final report - 12/15/75) 


b. North Dakota 


(l-’inal ri'i’ort - 5/li/7(i) 


4 - 


c. Other 


(Final report ~ 5/14/76) 


Determination of Kimcrical Const.nnts in SSI: Prediction Ht|uations 
includin{> Variable Significance Testing and SSIi Prediction 
Accuracy. 


Documentation for selected ITS for 


a. Kansas 


b. North Dakota 

c. Other 


(Final report - 12/15/75) 
(Final report - 5/l-l/”6) 
(Final report “ 5/14/76) 


Application of S.SE Prediction Equations to Cori'esponding Simulated 


Remote Analogue Areas for selected ITS 


Applications of tests for statistical differences between 
predicted and actual signatures in ficltls located in current 
ground dtstalcss or simulated ground dateless areas :.’hich arc 
identified by photo interpretation as being v.-hcat. Checks for 
classification accuracy v/hcrc ground data is actually available. 


1. Kansas 

2. North Dakota 

3. Other 


(12/15/75) 
( 5/M/76) 
( 5/14/76) 


b. Identification of optim.'’! SSH prediction inforir«at.lcn/plioto 
intcrj)rctcr remote who:.' field location feedback procedures. 

(PiU liminary assessment - 12/J5/7f>) 
final report - 5/14/76) 

Sensitivity Analysis in Selected Intensive 'lest Sites to Quantify 
SSE Prediction Performance Relative to Input T)qDCs, Quality and 
Qiiantity. (Preliminary assessment - 12/15/75 

final report - 5/M/76) 
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1.3 Task III; Variable Probability Saiuplinij for Acreage Estiwation (Task 4.1.3.1b, 

(KSRP, U.C. bcrkolcy) 


1.3.1 Introduction 

The human's ability to recognize wlieat growing areas and make quick 
estimates of wheat acreage using inexpensive pliotographic products can 
be incorporated into a variable probability sampling sclicmc. Uiider this 
procedure, sample units (segments) are selected with ]»robability proportional 
to estimated size (Pl’ES) of the acreage of wheat determiuod by ocular tcclmiqucs. 

1.3.2 Objective 

The objective of this task is to develop and test photo interpretation 
techniques for definition of v.hoat sampling strategy and to deCieo the PP1;S 
estimation formula associated v/ith the sampling strategy. The desired product 
is a variable probability sampling scheme utilizing im;igc interpretation that 
may be used, witliin specified application criteria, to allocate samples for 
unbiased wheat acreage estimation. 

1.3.3 Overall Approach 

The extent to wliich photo interpreters can be used to allocate siimples 
in an unbipsoJ and cost effective mariner for wheat .acreage cstiir.ation will 
be dotciTiiincd. To do this, seven subtasks will be accomplished. 

A. Photo Interpreter 1/iicat Estimates 

The ability of the photo intciq^rctcr to estimate wheat acreage 
of otlicr variables that can he correlated witli wiioat acreage and tlic 
extent of the corj*clnlion between the photo interpreter's estimate and 
computer estimate will be detcrm.i.ncd. lliis will be done in areas defined 
for R5RP LACin Study within t]ic U.S. that rcpi'escnt, as well as possible, 
the conditions encountered in the non-U. S. arc.'is. The ability of tiic photo 
interpretation scheme to handle the extremes in field size fouiid between 
specified areas such as found in the U.S.S.R., Cliina and India, will be 
specifically considered tlirough the use of U.S. analogues and the actual 
imagery from the countries. 

B. Definition o f PPES Estimation Formulae 

As a result of tlio investigation outlined in (A), detailed 
definition of the PPES estimators, tlic derivation of the cst ir.iviLors, tlic 
derivation of the estimators of variance, present at ion of sample size and 
sampling rate formulas, and the sample allocation luothoJ to moot the system's 
constraijits will be provided. 
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C. Estiiivition oC Snuplc Size 
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After determining the correlation between the cominitor,j>reclictions 
and human predictions of wheat acreage, the sample size requirements spec i fie 
to the v.'hcat growing regions under study by USUI* in the L\CII; will be estimated. 


D. Muman-Computcr Analysis of Segment Data 


If the human input proves cost-effective, a general pi*occdnrc for 
human and possible computer analysis of segments will be developed for the 
various wheat producing areas of the v/orld. This v/ill include the documentation 
of the advantages, as v;cll as the limitations, unusual situations that may be 
encoimtcrcd in the application of these techniques and tlic image types, and 
procedures required to estimate the selected i;hcat correlated parameters for 
spring, iifintcr, and mixed wheat areas. 


E. Effect of Hissing Data (minimum effort) 

Tlie best possible way to handle the problem of imagery which is 
missing and which is also important to tl)C PPES-plioto interpret at ion scheme 
will bo detenained. This procedure may entail the use of existing maps, 
meteorological data, and historical acreage data to generate surrogate ioagos 
to use in the photo interpretation estimation procedure. 

F. Change Detection (low priority) 

To the extent that any significant changes can 1)c detected within 
the cons traints of the existing se quent ia l data base for the intensive sites 
being studied, techniqiios will be evaluated to assess the impact of land use 
shifts on the sample allocation procedure. Tl\e PPES allocation using previous 
years of data will undoubtedly reduce tlie ability of the systcii to detect tiiis 
type of change. However, this problem can be overcome through a partial 
replacement sampling sclieme that will allow the inclusion of areas wlicrc cliangcs 
are likely to occur. 


1.3.4 Detailed Task Descriptions 


Tlie seven sub-tasks outlined in Section 1.3.3 will be conducted in a 
stepwise, coordinated manner. On-time ccmplctioii of all tasks will be 
dependent upon the timely receipt by the ilSRP of data from LACTE management. 
Tlie specific approaches for each of the seven sub- tasks will be as follows: 

^ • Photo Interpreter i '.'h oat Esti mates 


Optimal Inmian image-interpretation methods for wheat acreage 
estimation will he defined, described, and quantified. ’’Oplimnl" will 
be defined ns that coi;’binal ion of photo interpretation tcclinl([i!es and 
image sots allowing tlie liighcst correlations between ))!ioto inloj'prctcr 
wheat acreage estimates and "ku.own" (c.g. computer or ground truth) acreages 
within given cost ami time constraints. This research v/ill bo based on (1) 
the optimal imago intcvjirctation tcchniriuc comliinations defined in the KSRP's 
lAClE Tasks I (4,1.1.1a) and II (4.1.1.2f(l) and 4.1.1,2f(.3), and (2) investigations 
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to dctcriiriiic the optimal jji!a'»c sets alonp, v/itl\ photo iuierprctcr tcolmiqucs 
that will be coiiduclcd *,pecifically tor Task 111 (‘1.1.51b)* 

Optimal ijnanc-intcr]>retation technique combinations will be defined 
for major wheat producing types and spatial arrangement ))attcrns witliiii 
LACll: RSRP Study areas that rcprc.sont, as well a.s possible, conditions 
encountered in U.S. and non-U. S. areas. /\s a minimum, the ability of tlic PPliS 
photo interpretation scheme to handle extremes in field size and spatial 
arrangement bctv/cen tlic USSR, Cliina, and India will be ‘ipccifically considoroJ 
througli the use of U.S. analogues and the actual imagciy and supi)orting auxiliary 
data from the countries. However, in order to insure adequate evaluation of 
the photo interpretation techniques, tlic PPfS sclicmc, aiid llSUP's Task I and 
Task II techniques, data from all. countries involved in the L/.CIIJ CAperinent 
should be available for research. 


B . Definition of PPHS Estimation Pormulac 

Detailed definition of tlie PPIiS estimators, the derivation of tlic 
estimators, and the estimators of variance v/ill be conducted v/it!i a view to 
current and eventual LACIH system constraints. An example of a current 
constraint that would affect the matlicmatical representation of the estimators 
is the systematic nature of oligibic sample segment location im])OS'.;J by the 
GSFC system. Estimator bias and precision problems associated with such 
cin.T.*cnt constraints ’.vill be identified and contr.istcc! v/ith situations associated 
with alternative PPES derivations under less restrictive hACIE sampling segment 
locations guidelines. 

Derivation of sample size formulas will be presented jilong witli 
the overall sienplc design under wliich sa.mplc scgiiicnts may be allocated to 
specified I.ACII; countries, regions, zones, and acreage strata. Sai^plc size 
formulas and sampling designs will be specific to given LACIIi sy.slc],i coiw^traint 
sets. Countries, regions, and zones will bo as defined by I.ACU; ccoperators 
and management. Acreage strata will be defined on the basis of approximately 
constant levels of v/l\cat density as photo interpreted and checked with liistorical 
data so as to minimize the variatiori in wheat acreage estimates. 

C. Estimation of Sample Size 

For a given hACIE constraint set and associated PPES osSimators 
and scunplc size foriiiUlas, sam])lc segment sizes v/ill be calculated for Sjoccific 
country, region, zone, acreage strata within LACIF. IlSRP Study avca.s. This 
task involves acquisition (from JCS), indexing, preprocessing and data merging not 
already porrormed wnulor Tasks I and II for the areas o(; interest. In order to 
calculate sample sizes, acreage strata must be defined within I.ACTE reporting 
units, photo-interpreter estimates of whc.nt acreages must be obtained, and in 
some cases ground wlieat acreage data variance calculated wlicrc sucli i nformat Ion 
is available. In acMUion, Llio cost structure of the Fi’ES acreage estimation tcclmi(|ue 
will 1)C determined Including sami’lc unit costs at various stages or ph.asos in the 
sampling, acreage estimates will determine, in part, optimal sample r.itos for tlie 
rospocMve (intimation techniqin.is. Tlicsc calculations will ho reported for various 
regions at times sjiccified in R.SUP's I.ACIE Task d. 1,3.. lb milc.stono plan, Table d. 

Delay in 113RJ’ receipt of data will potentially delay sample siro riq)orting. 


RiTsoDOciBiLnY <2F ms 
ORfOINAL PAQB IS lOOR 


I) . 1 l unifin- Coi H]iutcr An;0 ysis of Sc :f>r,:cnt Data 

Based on tlic work in tho forego ivi« four sub-tasks, a f.oiicral 
procedure will bo outlined for coiiiliininy piiOl.u-intcr|irctaLion rosults and 
ccmputci* estimates in wljcat acrcaj'e ostimators. Tlio I’l’I-S r.ainplo dosijjn 
will ho specifically considered, thou'ih tlic application of Llie tocluiiqucs 
developed to date in other sample designs employed in LACIli will he 
considered as well. Included will be a documentation of advaiitagcs, limitations, 
and unusual situations to be encountered in tho use of the liumun-comimtcr 
segment .Dialysis for wheat acreage estimation. Image t>qics .and intei-prctation 
procedures required to estimate selected wheat correlated parameters for major 
v.iicnt producing area types considered in I-ACII; IISUP Study arc.ss v/ill also ho 
given. 

E. Effect of Missing Data (rainimiui effort) 


Wien missing HUTS imagery for given wheat production areas occurs, 
either due to cloud, hardv/arc, or processing problems, tccliniques will, be 
developed to generate cstiriU’tc-d photo-interpreter wheat acreage cstiii'atcs. 

-Tlic approach will be («) to generate and apply regression rclationsliips 
between previous plioto-interprcter estimates for the missing data ureas .and 
current plioto-intcrproter and/or computer estimates in simiJar areas where 
data is not currently missing; or (h) generate and ap))ly regression relation- 
ships betv/eon most recent photo- interpreter estimates and tlic following: 
existing map data, liistorical acreage data includiiig pr.:st photo- intciprcter 
and comjiuter estimates, current coinputcr estimates, and/or indlcu'S of malcoro- 
logical conditions. Apjilicalion of the regression relationship!-, ’./ill alio;/ 
least-square estimates of current photo-interpreter wheat acreage estimates 
in areas cun'cntly missing data. Co’iibincJ with this missing data approach will 
be a preliminary analysis of its effect on PFES e-Uimator cost, accuracy, 
and precision performance. 


F. Change Detection (low priority) 


It is presently conceived tliat coriplcto photo-interpreter 
estimates of wlic/it acreages v/ill be fortnulatcd only at tlic initial Invcntov-y 
period for those areas selected by LACIE inanagcmcat for subscipiont photo- 
interpreter Fl’f.S application. In order to inexpensively detect land use 
conversion to or from wheat production, it is projiusod that a plioto- interpreter 
sample clicck of v.iicat acreage estimates bo made at subycarly, yearly, or multiycarly 
intervals. Based on photo-interpreter results, sample segments for computer 
analy.sis can he reallocated liy a partial replacement strategy at .similar intervals 
in order to miniiiiize bias in the wheat acreage estiimitcs. 


Ail A c. a. v/ii j I 

change as it rclatc.s to 


UlllCiUIlL 


land u.'-.e 
xnmined. 


1.3.S Resource Data 
I. Existing mosaics 
II. lixi.sting mo.saics 
III. Existing mosaics 


Rc/juirements 

for State of Kan.sas 

for otlier states (Idaho or Montana?) 

for selected l-oj-oign Area (Oblast?) 


Date Renuived 
ay RSkP 

y/15/75 

10/1/7S 

in/ib/7.s 




Dat.c llw|uircU 
by iiSIU’ 


IV. . Updated state and country mosaics as they bocoi^'c available 
V. Selected UNDSAT color infrared composite .transparencies 

a) Kansas 

b) Other state (Idaho or Montana?) 

c) Selected Foreign Area (Oblast?) 

VI. Suppoi'ting cultural practice, climatic, wheat acreage, wheat 
variety, yield, soils, and other support data (T13D) for the 
above-named states and countries when not available through 
other LACIb U.C. BcrKolcy Tasks. 

a) Kansas 

b) Other state 

c) Selected Foreign Area 

VII. Field identification data, including SUS and other available 
information for above-named states and countries when not 
available thrcv.gli otlier L\Cin U.C. Berkeley Tasks. 

a) Kansas 

b) Other state 

' c) Selected Foreign A.rca 

VIII. Detailed description of capabilities and rcstrictiojis of 
hardv.'ure and software components (existing and developing) 
available to lACIF and JSC and Goddard 

IX. Complete description of the cm'rent LACIH sajjipllng plans 
for the USA and other LACIl; countries 

X. Description of other sai;ipling U fi D approaches, locations 

of application (c.g. particular ai'cas .and reporting district 
types of Kansas), and estimator results and variances. 

XI. Crop/land use (wheat and other crops/lnnd use types v:hci*c 
available) jn'ojjorv. i ons from computer class i.fic;ition for 
segments in areas defined for study. Also, co;.iputor eJassifie- 
ation accuracies by segment by crop/Iaml use type where availa! 

a) Kansas 

b) OLlicr state 

c) Selected Foreign Area 



as .available' 


7/ir./75 

10/1/75 

lO/l.S/75 


7/15/75 

10/1/75 

10/15/75 


7/15/75 

lO/J/75 

10/15/75 


7/15/75 


7/15/75 


7/15/75 


)Ic. 

7/15/75 
10 / 1/ 75 
10/15/75 


— . 1 


TAni.i: A MTLnr;ro;>ii: plan i-ou lacip. task ni. 

1, DctcrniJnation of Photo Intoi))rotor V/lioat. list j mat ion Abilities for 
selected Worth Ajucrican (ih\) V/Iioat Prodiieinj; Areas. 

a. Kansas (counties) (Kcsultr. 12/S1/7S) 

b. Idaho (North Uahota?), Other state (fi/1‘1/76) 

c. Selected Porcign Area (5/l'l/7G) 

2, Definition of PI’S KstiHation Pormulac 

Includes the sample design for the wheat acreage ostiipation 
problem along \.’ith the associated estimators, estimators of 
variance, sajiiple si"c and rate formulas, and the specification 
of wheat acreage stratification strategy. 

(Preliminary documentation - enough to iisplone/Jt 
plan - 12/;il/75 

Final report - .S/ld/76) 

3, Fstimation of Sample Size 

Includes any additional data acquisition from JSf' and data ju'o- 
ccssing not already specified by Tasks 1 and II; also includos 
acreage stratification, plioto interpreter acreage ost Jinnt ion, 
selected ground data varia-iice analysis, sampling cost analysis. 


/jialysis for selected Nortli American areas; 

a. Kansas 

b. Other states (Mortli D.ikota, Idaho??) 

c. Oblast? 


(12/31/75) 
( VM/7G) 
( 5/14/76) 


4 . Ucsults of V.lieat y\creage T-stimaticn Including proportions and 
variance of prtjportions and related str.ti sties suiriMarized by 
counties and/or other reporting districts tliouglit to be 
desirable. Included also, if necessary, arc jiiri)s of counties 
and/or other reporting districts showing acreage str.’.ta and 
associated proportion estimates. 


a. Kansas 

b. Other state 

c. Oblast . 


(12/31/75) 
( 5/14/7G) 
( 5/14//G) 


5. Develop a General Procedure aiid Criteria for Human and Coinputer 

Analysis of Sanplc Segments (5/14/76) 

6. Del(Mininc Mothed:; to Develop Plioto Interpreter hsi i mates i'dien 

Missing Data Occurs (5/14/76) 

7 . Dcvoloi) tecliniqncs to Detect Land Use Ch.ange Affecting Wheat 

Acreage Hstimales . (5/14/76) 









U.C. BERKELEY SIGNATURE EXTENSION APPROACH 

(1) Stratify regions into analogue areas tending to give rise to 
a similar sequence of wheat signatures. 

(2) Identify land cover signatures in recognition segments by 
utilizing land cover, phenological phase specific signature 
prediction regression equations developed in training segments. 

(3) Classify recognition segment based on training signatures from 
the recognition segment identified by the process in (2). 
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STRATIFICATION 


Analogue areas are strata of specified extent defined by 
relatively homogeneous combinations of values for environment 
variables such as 

® General Soil Type 
®Land Use/Cropping Practice 
® Growing Season Day-Degree Sum 
® Growing Season Precipitation Sum 
® Average Last L-ate of Spring 
Frost 

® Average Temperature for the 
Coldest Month 

Analogue area boundaries may be updated when desired, generally 
every several years. 




i 
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Spectral Signature Statistics for General Cover 
lypes Within Recognition Segments 

(1) Within a given analogue area type, signatures are extended between 
the training segment (s) and the recognition segment (s) by use of 
regression equations. 

(2) There is one regression equation developed from training segments 
for each combination of major cover type, major phenological phase, 
and LANDSAT channel. 

(3) Each regression equation is designed to predict the spectral signature 
in a given LANDSAT band for a given cover type according to predictor 
variables derived from the recognition segment such as 

o Crop Maturity Codes from A.I.s 
o Soil Color Code from A.I.s 
o Growing Season Day-Degree Sum 
0 Planting Season Day-Degree Sum 
o Growing Season Precipitation Sum 
® Planting Season Precipitation Sum 
o Sun Angle 
o Scan Angle 

o In-Segment LANDSAT Spectral Data Such 

as the 7/5 Ratio (Crop Calendar Related) 

» Atmospheric Transmittance (Generated from 
Ground Horizontal Visibility Data) 

(4) Predicted signatures are statistically matched to fields or clusters 
in the recognition segment 

(5) The spectral statistics from the matched fields or clusters in the 
recognition segment are then used for recognition segment discriminant 


3 / 


analysis. 
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SPACE SC1I.KCC5 LAliOllATOnY 


i»;mu:Lr.Y, camkoilvia 94720 


•Tl'/X: UC SPACE BERK 
(910) 366-79^S 

Remote Sensing Resenreh Program 
260 Space Sciences Uboratoiy 
University of California 
Berkeley, California 94720 

Hr. Joe Garcia July 8, 1975 

Code TF3 

KASA“Jolinson Spacecraft Center 
Houston, Texas 77058 


Dear Joe: 

Enclosed please find a list of the U.C. Berkeley P.cnotc Sensing 
Research Program's (P.SI'P) research and development data needs for the kAClE 
spectral signature extension task. It would be desirable to have as r.ucli of 
this data as possible for tlie LACIE segment sites as well as fur tlie inten- 
sive test sites. It is expected that this data will be provided througli 
task II when not available or duplicated by the other RSRP I.AC1E task's data 
requests or when not readily available through public sources. 

Some of the items on tlie list would theoretically not be used in an 
operational system, but arc necessary for the research and development of 
the system. 

Some of the data items arc already in l^and as you know, but are in- 
cluded for a full documentation of the data requirements for this task. 

Also I have rnde notes on the availability from public sources of certain 
items of data so that if other researchers have similar data needs you can 
advise them of their ready availability through public sources. 
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REPRODUCBILITY OF THE 
ORIGINAL FAOC D POOR 


KSUP lACIl; Sij’.nnlu rc nxlc ns ion K fj 1) Data N eeds* 


1. County liisioricnl crop data 

Partly nvail.nhlc in U.S.D.A. Crop and Livestock Reporting; Service's 
annual crop acrcaijc summaries by counties. Need yct^rs '71 through 
•75 for all tost sites. 

2 . Crop calendar information 

Normal crop calendars for the test site areas with adjiistncnts for 
the '73- '74 ^ '74 -'75 seasons. 

3. Cultural practice infonaation for specific regions witliin the study: 

a) amount of irrigation practiced, v;hat system used, location map 
of irrigation projects within the regions of the test sites. 

b) crop rotation practices, fallow periods, double cropping or 
overcropping practices specific to given crop types and .ns 
specific to gcograpliic regions such as crop reporting disTricts 
as possible. 

c) type and prevalence of fertilizer applications by specific 

crop or soil conditions. Similarly for insecticide and herb- 
isidc treatments. , » 

d) names and characteristics of major varieties of specific crops 
as specific to geographic regions (crop reporting districts) as 

• possible. 

e) significant variation in seeding rates and plant densities due to 
different planting practices by crop as specific to region as 
possible. 

4. Base map photography of test sites 

5. Any e.xisting LAIxDSAT coincident aircraft pbotograpliy with coincident 
or nearly coincident ground data for test sites. 

6. IJ^NIJSAT color infrared, 9" X 9" transparencies to cover the entire 
state (s) of study. For Kansas one each within the foilwing time 
periods: 

a) 1st week of July 

b) 1st week of August 
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7 . 

8 . 


c) i.'iiil-Octol>cr (;im1/or mul-March) 

Sclcclcil CCTs to cover test site locntions 
Soils data coiisistinc of; 

a) the 7th ;i}3i>roxi nation classification to the subnroup level by 

state. 

b) Local soils maps to the series level, (these arc j'artly available 
in the U.S.n.A. S.C.S. soil survey reports for specific count ic.s, 
however, not all counties have been mapped as yet or tlie reportr- 
arc at this time not easily available to the j)ublic). 'fliesc 
reports shonlil contain information on the profile characteristic:; 
of the soils, the available water holding capacity of sjiccific 
soils, etc. 

9. Climatic data as specific to test site locations as possible. Tne 
data should include daily tem]tcratures, daily i.riximum tci^iper.'Uu.-c, 
daily ninimuia temperatures, daily precipitation, and d.-’.ily Iln (all 
wave net radiation input), lliis data is partly available tbrnugli 
the Department of Conferee, r.nvircninental Data Service’s (Aslievi 3 le, 
N.C.) monthly state cl imatalogical data reportr.. 

10 . Atmospheric data as specific to test site locations ns jjossible. 

Tne data should include: • 

a) horizontal visibility data such as reported by aij'port weather 
observation stations. 

. b) sky cover * 

c) *i possible sunshine 

d) wind vector 

11 . Ground data for test sites specific to individual fields consisting of: 

V- 

a) crop identification 

b) variety 

c) stand (luality with specifics of any stress conditions (cau;;cs of 
stress) 

12. information--any publications dealing t.'ith crop bobavior 
and iminagemcnt requirements by variety. 


Data requested should be provided when not .available llirougb other Kbl'.P 
LAClil tasks or individu.il KSRP effort. 



* 


